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Abstract 

We discuss some properties of a supersymmetric matrix model that is the dimensional 
reduction of supersymmetric Yang-Mills theory in ten dimensions and which has been 
recently argued to represent the short-distance structure of M theory in the infinite 
momentum frame. We describe a reduced version of the matrix quantum mechanics 
and derive the Nicolai map of the simplified supersymmetric matrix model. We use 
this to argue that there are no phase transitions in the large- limit, and hence that 
S-duality is preserved in the full eleven dimensional theory. 



Work supported in part by the Natural Sciences and Engineering Research Council of Canada. 



The conventional understanding of the spacetime structure of string theory has drastically 
changed over the last few years. It has been realized recently that all ten dimensional su- 
perstring theories are related by non-perturbative dualities and that they can be thought of 
as originating, via Kaluza-Klein types of compactifications, from a single, eleven dimensional 
theory known as 'M theory' (see for recent reviews). The dynamics of this theory are not 
yet fully understood. Some of the central objects in understanding the string dualities are 
non-perturbative, p-dimensional degrees of freedom known as Dp-branes , on which the end- 
points of strings can attach (with Dirichlet boundary conditions). The low-energy dynamics 
of a system of parallel D-branes can be described by an x matrix model obtained from 
the dimensional reduction of ten dimensional supersymmetric Yang- Mills theory with gauge 
group U{N) 0, The large- limit of this matrix model has been recently conjectured to 
describe the small distance spacetime structure of M theory in the infinite momentum frame 
P]. The explicit solution of this matrix model therefore has the potential of providing a 
non-perturbative description of the largely unknown dynamical objects describing the short 
distance behaviour of the full eleven dimensional theory 0. 

In this Letter we will discuss some basic properties of the supersymmetric matrix model 
introduced in |^. We examine a particular reduction of the model proposed in |^ 
to static D-brane configurations with transverse 5*0(8) rotational symmetry. We explicitly 
construct the Nicolai map associated with the supersymmetry in this reduced model and use 
it to analyse Schwinger-Dyson equations of the matrix model. We show that the results from 
this analysis are consistent with other known results of the D-brane field theory, and also 
that this simple approach gives some insights into the structure of the full matrix model. In 
particular, the reduced model seems to have no phase transitions in the large-A^ limit and 
S-duality is preserved in this representation of the full eleven- dimensional theory. 

First, we discuss some aspects of the representation of systems of D-branes by Yang- Mills 
fields. Consider the gauged supersymmetric matrix quantum mechanics with action 0, || 
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where DtY = dtY — z[Ao(t), Y] is the temporal component of the gauge covariant derivative, 
the trace is taken over unitary group indices, and we have chosen units in which the string 
tension is a' = l/2n. Here X'{t) = [X^fe(t)], a,b = 1,...,A^, i = 1,...,9, are A^ x A^ 
Hermitian matrices in the adjoint representation of U{N) which are obtained as the spatial 
components in the reduction to (0 + 1) dimensions of a (9 + 1) dimensional U{N) Yang-Mills 
field A^{x, t), ^ = 0, . . . ,9. They describe the collective coordinates of a system of A^ parallel 
DO-branes (with infinitesimal separation), and they transform under the vector representation 
of the rotation group SO (9) of the space transverse to the compactified eleventh dimension of 
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the underlying supergravity theory. The superpartners of the matrices X* are the Majorana 
spinors ^/'"(t) = ["^afel^)]' ^ ~ 1, • • • , 16, which transform under the 16-dimensional spinor 
representation of S'0(9), and under the adjoint representation of the gauge group U{N). The 
Dirac matrices 7j are the generators of the spin(9) Chfford algebra 

{7.,7,} = 25., (2) 

in a Majorana basis. The coupling constant g is related to the eleven dimensional compacti- 
fication radius R by 

R = 9"'lp (3) 

where Ip is the eleven dimensional Planck length. 

The action (|I]) describes the short-distance properties of DO-branes in weakly-coupled type- 
IIA superstring theory 0, It was argued in [^] to be the most general infinite momentum 
frame action with at most two derivatives which is invariant under the U{N) gauge group 
and the full eleven dimensional Lorentz group. It is further invariant under the infinitesimal 
Af = 1 supersymmetry transformation 

= M E Dtx\^r, + 1 E ^^1 7,]? ) (4) 



5eAo = -2£>„ 



where are 16 global Majorana spinor parameters. There are another set of 16 super- 
symmetries which are realized trivially as Se'ip"' = ^e'X^ = ^e'^o = 0, where e and e' are 
independent supersymmetry parameters. Together, these two sets of supersymmetry transfor- 
mations (with slight modifications) generate the full 32-dimensional 16 ©16 representation of 
the super-Galilean group of the eleven dimensional theory in the light-cone frame. A number 
of properties of M theory have been verified using the action (|I]) 

It may seem puzzling at first sight that a dimensionally-reduced gauge field A yields the 
appropriate D-brane coordinatization. What is intriguing though is that string theoretic T- 
duality is a key to this connection. If we compactify the first spatial dimension onto a circle 5*-^ 
of radius TZ (or more generally several dimensions onto a torus), then the angular coordinate 

(or coordinates x*) takes values in the interval G [0, 2ti7V\. The original supersymmetric 
Yang- Mills theory in ten dimensions describes the low-energy sector of open superstrings, and 
in this theory there exists other gauge-invariant variables, namely the path-ordered Wilson 
loop operators 

W[A] = tr Pexp Ai dx'^ (5) 
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which are invariant under the large gauge transformations which wind around the compactified 
direction. The connection appearing in the argument of the exponential in (^) lies in the 
adjoint representation of U{N). The corresponding classical gauge field orbits are topological 
and lie in the interval Ai G [0, For each coordinate living on a circle 5*-^. there is a "dual" 
coordinate X* = a'Ai which also lives on a circle S^^ but with a dual radius = a'/TZi. As 
we have mentioned the ten-dimensional supersymmetric Yang-Mills theory describes the low- 
energy dynamics of open strings with Neumann boundary conditions. Under T-duality TZ —>■ 
a'/TZ the Neumann boundary conditions are transcribed into Dirichlet boundary conditions. 
But this means that the topological gluon degrees of freedom are converted into the D-brane 
fields X{t) = a'A{t) describing open strings with Dirichlet boundary conditions. Thus the 
ten-dimensional gluon fields describe the dual theory of the D-branes and T-duality naturally 
identifies the topological orbits of the Yang-Mills fields as the D-brane coordinates. The role 
of T-duality in the matrix model (|l|) has also been addressed from other points of view in p|. 

A question which now arises is how to really measure the original angular coordinate 
given the gauge field coordinate X\ The solution is to consider more complicated objects, 
such as the abelian Wilson loops with a winding number n G Z 




which are associated with a single D-brane configuration. In the general case, i.e. in the case 
of a non-abelian gauge group, we should consider Wilson loops in different representations of 
U{N), but for the sake of illustration we shall discuss here only the U{1) Wilson loops (^. 
This simplication can be thought of as a compactification of all of the ten dimensions on which 
the group variables are restricted to the maximal torus of the U{N) gauge group. Using 
these objects and T-duality one can in principle obtain the angular coordinate through 
a superposition J2n£z'^n[A] of different Wilson loops (using harmonic analysis on unitary 
groups in the general case). 

In fact, this construction demonstrates that the Wilson loops in the "N-theory" (i.e. ordi- 
nary open strings with Neumann boundary conditions) are equivalent to the vertex operators 
in the "D-theory" (open strings with Dirichlet boundary conditions). To see this, we ex- 
pand the topological gauge field configurations of the low-energy description of strings on the 
compactified space in Fourier modes as 

A,{t,x)=Y: 4'^ (t) e^^^^ , xe[0, 2n/n] (7) 

These dual fields are described by a (1 + 1) dimensional gauge theory that is also dimensionally 
reduced from the ten-dimensional supersymmetric Yang-Mills theory |^. Then the Wilson 
loop (^ becomes 
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where Pn = n/r = nTZ/a' is the momentum of the string winding mode in the compactified 
direction. Thus T-duahty converts the non-commuting position matrices X* describing the 
D-brane configurations into gauge fields in the dual theory, and also the winding number of 
the topological gauge field modes into the string momentum in the compactified direction. 
In the general case then, we can conjecture the equivalence between the string scattering 
amplitudes of the D-theory defined by correlators of the vertex operators and expectation 
values of the Wilson loop operators in the dual N-theory0, 



{W^,[A]---W^^[A]) 



N 



D 



(9) 



where fij G and we have defined Vr,^^. [A] = ni*=i W^(raj)»[^] terms of the abelian Wilson 
loops (P). In the non-abelian case, the correlators @j will generalize in the appropriate way 
in terms of representations of U{N). 

We now examine the problem of obtaining an explicit solution of the matrix model 
For this, we further dimensionally reduce the theory described by (|^) to a zero- dimensional 
N X N super symmetric matrix model, i.e. we ignore the time dependence in (|I|) and work 
in the Weyl gauge Aq = 0. This means that we are studying the model separately over each 
constant time slice describing a static configuration of the DO-branes. This reduction can be 
thought of as originating by compactifying the time direction of the ten dimensional Yang- 
Mills theory where the adjoint representation fermions have periodic boundary conditions, and 
then taking the limit in which the radius of compactification vanishes. This simplification has 
the advantage of eliminating non-local operators that would appear from the time-dependence. 
We shall discuss the inclusion of time-dependent fields at the end of this Letter. 

With this further reduction, the partition function of the model is given by the finite- 
dimensional matrix integral 



nn^^:^ n« ^xp^-tr 



(10) 



We can expand the matrix integration variables in (|TOD in a basis of the unitary group 
as X* = X\T^ and ip"" = ip'^T^, where A = 1, . . . , N"^ and the Hermitian U{N) generators 
satisfy 
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AB rpC 



tr T^T 
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The integration over the Majorana fermions in (^) is Gaussian and can be evaluated explicitly 
using the Berezin integration rules for t/^^. It produces a square root of the determinant 
determined by the representation of the adjoint action of X\ and ([lOD becomes 

N 



cn 



iV2 

UlldXh 

i D=l 



Pfaff 



2g 



/^^^E(7.)^^, 



c 



(12) 



^See for another discussion of the relationship between string vertex operators and Wilson lines. 
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where the Pfaffian is taken over both the adjoint U (N) representation indices A, B = 1, . . . , 
and the spin(9) indices a, /3 = 1, . . . , 16. Here and in the following cn denotes an irrelevant 
numerical constant. 

We would now like to exploit the supersymmetry of the zero dimensional model to com- 
pute correlation functions of the matrix model. When the number of bosonic and fermionic 
degrees of freedom are the same, the supersymmetry is maximal, in that it holds even when 
the fields are off-shell. In the present model, the number of fermionic and bosonic degrees 
of freedom do not match. Normally, supersymmetry would require on-shell fields and the 
addition of auxilliary fields to make the number of physical boson and fermion modes equal. 
However, we can adjust things to match by exploiting the original interpretation of the matri- 
ces from the dimensional reduction of ten-dimensional supersymmetric Yang-Mills theory. 
The latter theory can be gauge-fixed and quantized in the light-cone gauge, after which there 
are only eight propagating gluon degrees of freedom corresponding to the various possible 
transverse polarizations. Since a Majorana-Weyl spinor in ten dimensions has eight physical 
modes, the minimal Yang-Mills action in ten dimensions is supersymmetric without the need 
of introducing auxilliary fields. Thus we match the collective degrees of freedom X* of the sys- 
tem of D-branes with the physical modes of the full Yang-Mills theory by reducing the target 
space degrees of freedom from nine to eight by setting = and working in the Majorana- 
Weyl representation of spin(9) in the matrix model above. The constraint X^ = can be 
thought of as a light-cone gauge fixing condition in the nine-dimensional transverse space. 
Although not precise from the point of view of the M theory dynamics, this simplification 
produces a toy model that will shed light on some of the properties of the nine-dimensional 
theory that we started with[|. 

With this simplification we now exploit some features of the group theory for S0{9). The 
Dirac generators of spin(9) in the Majorana-Weyl basis can be constructed from the reducible 
8s © 8c chiral representation of spin(8) by decomposing the 16-dimensional gamma-matrices 
in the 8x8 block form 



•^Some different reductions of the matrix quantum mechanics (|^) have also been suggested. In jT^ it was 
argued that the model can be truncated to zero dimensions by augmenting the transverse rotational symmetry 
to 50(11) and viewing the matrix model as the dimensional reduction of supersymmetric Yang-Mills theory in 
(10 + 2) dimensions. In jlOj it was argued that the D-brane field theory associated with weakly-coupled type- 
IIB superstrings could be viewed as the large- A'^ reduction of the ten dimensional supersymmetric Yang-Mills 

theory. 
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and similarly with dotted and undotted cliiral indices interchanged. The spin(8) Dirac gen- 
erators can be expressed explicitly as direct products of 2 x 2 block matrices by 



(71 )a = -^Cr2 (g) 0-2 ® (T2 
(73)" = il(g) 0-3 (g) (T2 
(75)0 = i(T3 (g) 0-2 ® 1 
(77)" = 



(15) 



= il (g Oi ® 0-2 
= iai (g 0-2 ® 1 
= 202 ® 1 (8 0-1 
= 1 ® 1 ® 1 

are the Pauli spin matrices. The 
(1 

remaining spin(9) Dirac matrix is then 79 = 7172 ■ ■ ■ 78 = 

Vo -1, 

The block decomposition (0) of the first eight gamma- matrices shows that the Pfaffian 
in ( [T^ ) with = becomes squared, with the spinor part of the determinant restricted to 
the spin(8) chiral indices. Then we can write the partition function as 



where cxi 




and 03 



cn 



8 iV2 

UUdXh 

i=l D=l 



det 

A,_B;l<a,d<8 



1 



^/^^^E(7.)^^c 



exp 



'^9 l<i<j<8 



(16) 



This effective reduction to eight dimensions matches bosonic and fermionic degrees of free- 
dom and allows us to exploit the supersymmetry in a simple way to completely solve the 
matrix model. Essentially it enables us to use the triality property of the eight-dimensional 
rotation group, i.e. that there exists automorphisms between the vector and chiral spinor 
representations of 5*0(8). 

We now label the eight spatial indices i as the chiral indices a, d of the spinor represen- 
tation, and the first eight components of the 16-component spinor field ip as the chiral parts 
il)"" and the last eight components as the anti-chiral parts in the 8s © 8c decomposition 
of spin(9) above. The static reduction of the quantum mechanical action (P can then be 
written in the standard form of an TV = 1 supersymmetric field theory (after integration over 
superspace coordinates) as 



2^ 



E 

i,A 



dF 
dxi 



dXldX%^ 



where the super-potential is 



= 3(7.)° tr X^ [X",X. 
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d^F 
dX%dX%^ 



B 



(17) 



The representation (0) can be derived using the symmetry properties of the gamma-matrices 
([l|), the U{N) Jacobi identity 

yABC j!ADE _ j^ADB jACE _ j^DC jABE j^j^g^ 
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and the 5*0(8) Fierz identity 



(7.)^(7,)^ = M« + (7.,)^ 



(20) 



where we have introduced the spinor matrix 



a 



) 



(21) 



and similarly for {'Jij)^- 

The form of the action (|l^) identifies the Nicolai map (i.e. the Hubbard- Stratonvich 
transformation for the bosonic potential J2i<j^''^[^^j^''Y) IHI of this supersymmetric field 
theory as 



that comes from integrating out the chiral fermion fields. The partition function is thus 
trivially a Gaussian Hermitian matrix integral and is formally unity, 



The free energy logZ is thus trivially an analytic function of the coupling constant g every- 
where and it does not exhibit any phase transitions, even in the large- limit. Furthermore, 
the correlation functions which are invariant under the supersymmetry transformations (^ 
can be obtained by differentiating the free energy with respect to the coupling constants of 
the model (in an appropriate superspace formulation). Thus any supersymmetric correlator 
of the model vanishes, which is just the standard non-renormalization that usually occurs in 
supersymmetric field theories. The existence of the Nicolai map and these implied properties 
of the matrix model are essentially the content of the supersymmetric Ward identities. 

The only observables of the matrix model which are non-trivial are those which are not 
supersymmetric. To examine such correlation functions, we use the Nicolai map (^2]) to ex- 
press correlators (■) of the original matrix model (normalized so that Z = 1) as free Gaussian 
averages {{■ )) of the Nicolai field. For instance, from ((W^"^)) = we deduce {[X^X^"^) = 0. 
This means that the classical ground state of the model (the minimum of the bosonic potential 
X]i<j tr[X*, X-']^) is that wherein the D-brane coordinates commute and have simultaneous 
eigenvalues corresponding to definite DO-brane positions. The full matrix model, which in- 
corporates quantum fluctuations about the classical ground state, thus describes smeared-out 
DO-brane configurations in a spacetime with a non-commutative geometry |^, |^ . 



W,^{X) = ^ = 2 (7.)«r'^c^^Xf or W^f = -WX^, X^r (22) 



From ([TTp we see that the Jacobian factor | det[9PF^/(9X;g]| ^ which arises in the change of 
variables X W{X) in the partition function (|16]) will cancel exactly with the determinant 




(23) 
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More generally, we note that the Nicolai map X W{X) is many-to-one, so that general 
correlators of the X matrices can have a multi-valued branch cut structure. To see if this is 
the case, we use the Nicolai field to write down a set of Schwinger-Dyson equations for the 
matrix model. The basic identity follows from the formula for Gaussian averages of products 



of even numbers of the fields VTf'', 



TI/a2n-l62n-l TI/a2„62n 
i'2n-l i2n 



3_ 

N. 



[Saib2^a2biSa3b4Sa4bi " " " Sa2„b2n-l + n^"'^-* [oi, bi] 02, 1)2'-, ■ ■ ■ '■, a2n, &2 



(24) 

where 11 contains the sum of delta-functions over all permutations of indices. The delta- 
functions in the indices ik come from the 50(8) invariance of the reduced matrix model, 
while those in the indices a^, hk arise from U{N) invariance. The non-vanishing correlation 
functions of the model are those which respect both of these symmetries. 

As an explicit example, we set = 2 in (53) and sum over ii = 12, is = ^4 and ai = 625 

a2 = bi, 03 = 64, a4 = 63 to get 



tr 

N 



{ 1 



2 

iV2 



(25) 



In the large-X limit, the expectation value of a product of invariant operators factorizes into 
a product of correlators. Thus at X = 00 (pSf ) becomes 

2 



E 



/tr 

\n 



X\X^ 



10^4 



(26) 



On the other hand, setting ii = ^3, ^2 = H and the a's and 6's equal in the same way as above, 
we get the N = 00 equation 



x\x^ 



(27) 



Combining (126|) and (^) together we find that the large-X invariant variance of the 5*0(8) 
operator X]^ is trivial. 



tr 



x\x^ 



x\x^ 







(2^ 



Note that (|28|) is a stronger statement than just the large-X factorization of correlators, as it 
implies a non-trivial factorization over the 5*0(8) indices as well. This is one manifestation 
of the supersymmetry of this matrix model. Similar other such identities can be derived for 
higher-order correlators of the X-fields. In this formalism, it is also possible to treat n-point 
connected correlation functions of the model. 
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The forms of the correlators above (and in particular (^)) seem to suggest that all 
non-vanishing observables in this model are analytic functions of the coupling constant g 
aX N = oo. This in turn implies that the large- limit of the matrix model exhibits no phase 
transitions as one continuously varies g. As the relationship with M theory dynamics is even- 
tually obtained in the uncompactified limit where i? ^ oo the Nicolai map demonstrates 
explicitly that this limit can be taken unambiguously since there is no variation in the ana- 
lytic structure of the large- solution. Moreover, the absence of phase transitions suggests 
that S-duality g 1/g is maintained in the large- A^ limit of the matrix model above. A 
more precise examination of these properties requires the inversion of the Nicolai map (^2]) 
to get X{W), which would enable one to compute arbitrary non-supersymmetric correlators 
of the matrix model. The problem in trying to construct this inverse map is that generally 

7^ 0, corresponding to the fact that the D-brane coordinates live in a non-commutative 
spacetime, so that it is not possible to simultaneously diagonalize the X*'s and find the rela- 
tionship between the eigenvalue models for the W and X fields. The entire non-triviality of 
the matrix model lies in the correlators of invariant combinations of the operator X{W). The 
problem of inverting the Nicolai map has been discussed from a perturbative point of view in 
p!3| , where it was also shown that this transformation is a non-polynomial functional of the 
bosonic fields in the ten dimensional M = 1 supersymmetric Yang-Mills theory. It would be 
interesting to determine this inverse map, and use it to examine the properties of the Wilson 
loop correlators that we discussed earlier in addition to the large-A^ analyticity features of 
general correlators of the matrix model. In any case, we have formally described a method in 
which one can study features of the string scattering amplitudes (|]). 

The results described above are only precisely valid with both the elimination of the tem- 
poral dimension and the reduction to an 5*0(8) spacetime symmetry group. If we reintroduce 
the time dependence of the matrix variables then the Nicolai map is determined as the non- 
local, time-dependent functional W{t) = DfX'{t) - W'{X{t)) with W' given in (0). Then 
the partition function yields the winding number of the multi- valued Nicolai map. The above 
results from the reduced matrix model, such as the analyticity in the coupling constant g, and 
hence the S-duality in the eleven dimensional compactification, show that the simple method 
discussed above has the potential of providing some insights into the structure of M theory. 
It would be interesting to see if the reduced matrix model can describe other features, such 
as membrane interactions 0, of the eleven dimensional theory. It would also be interesting 
to determine if the Nicolai map obtained above can be used to describe any properties of the 
ten dimensional supersymmetric Yang-Mills theory itself. 

We thank O. Lechtenfeld for comments on the manuscript and N. Mavromatos for inter- 
esting discussions. 
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